Estimating herbaceous aboveground biomass in Sahelian rangelands using Structure from Motion data collected on the ground and by UAV

Abstract Herbaceous aboveground biomass (HAB) is a key indicator of grassland vegetation and indirect estimation tools, such as remote sensing imagery, increase the potential for covering larger areas in a timely and cost‐efficient way. Structure from Motion (SfM) is an image analysis process that can create a variety of 3D spatial models as well as 2D orthomosaics from a set of images. Computed from Unmanned Aerial Vehicle (UAV) and ground camera measurements, the SfM potential to estimate the herbaceous aboveground biomass in Sahelian rangelands was tested in this study. Both UAV and ground camera recordings were used at three different scales: temporal, landscape, and national (across Senegal). All images were processed using PIX4D software (photogrammetry software) and were used to extract vegetation indices and heights. A random forest algorithm was used to estimate the HAB and the average estimation errors were around 150 g m−² for fresh mass (20% relative error) and 60 g m−² for dry mass (around 25% error). A comparison between different datasets revealed that the estimates based on camera data were slightly more accurate than those from UAV data. It was also found that combining datasets across scales for the same type of tool (UAV or camera) could be a useful option for monitoring HAB in Sahelian rangelands or in other grassy ecosystems.


| INTRODUC TI ON
Grass or herbaceous vegetation is present in most ecosystems on earth, and herbaceous aboveground biomass (HAB) has been used as an indicator of ecosystem productivity and functioning in a plethora of studies in relation to biodiversity-ecosystem functioning relationships (Hector et al., 1999). HAB is also the main source of feed for many wild and domestic animals, and is used in many studies to assess feed availability Ickowicz, 1995).
HAB is usually collected by destructive sampling (removal of vegetation), which is a laborious work method, so the development of indirect sampling using remote sensing methods has been under continuous development for several decades (Reinermann et al., 2020). Many tools are based on measuring the spectral characteristics of vegetation. Such measurements are normally used to calculate indices (proxies of vegetation) based on a combination of near-infrared and visible reflectance of the vegetation (Bannari et al., 1995;Barbosa et al., 2019;Rouse et al., 1974). These indices are generally based on satellite remote sensing images, but can also be derived from images directly taken on the ground with a camera, or from near-field remote sensing using Unmanned Aerial Vehicles (UAV) (Candiago et al., 2015;Cruzan et al., 2016;Grenzdörffer et al., 2008). Another tool for estimating HAB using remote sensing is the LiDAR technique (Schulze-Brüninghoff et al., 2019). Therefore, the volume of herbaceous vegetation is subsequently derived from a LiDAR-retrieved point cloud, which can be used as a proxy of HAB.
Some studies have tested combining both the volume of the herbaceous vegetation and the spectral signature to estimate HAB .

Structure from Motion (SfM) is a photogrammetry process that
can be used to create a variety of 3D spatial models and 2D orthomosaics from images (Schonberger & Frahm, 2016). The process relies on a set of images of the same object taken from different angular views, and the SfM produces a 3D point cloud of the object. This process is generally applied to UAV images and the point cloud from the set of images is used to produce an orthomosaic and a digital surface model (Frey et al., 2018). Color indices can be calculated based on the orthomosaic depending on the sensors used (multispectral or simple RGB sensors). Applying UAV imagery with the SfM process is widely used to estimate HAB from 3D-based indices and vegetation indices (Aasen et al., 2015;Cunliffe et al., 2016;Kolarik et al., 2020;Possoch et al., 2016;Wijesingha et al., 2020). The SfM process can also be applied to imagery generated directly from sensors mounted on the ground. This process has been used to develop a 3D model of individual plants under greenhouse conditions, or directly in the field Andújar et al., 2018;Bossoukpe et al., 2020;Cooper et al., 2017).
Most studies using the SfM process have been carried out in temperate grasslands and on a local scale (Lussem et al., 2018(Lussem et al., , 2019Wijesingha et al., 2020). In this study, we assessed the applicability of the SfM methodology for estimating HAB in Senegalese savannah ecosystems including a large diversity of plant functional types. The objectives were to: • Assess the applicability of the SfM methodology based on reflectance in the visible part of the spectrum (red-green-blue; RGB) for biomass monitoring.
• Compare the biomass estimated using the SfM methodology based on images collected both with UAVs and with a digital camera on the ground.
• Study whether the SfM methodology can be used to monitor HAB variability over the growing season, across a savannah landscape locally, and spatially across a diversity of plant functional types.
• Lastly, compare calibration carried out across the different datasets and calibration within each dataset in order to see whether existing datasets can be used for such calibration.

| MATERIAL AND ME THODS
We collected RGB images using a camera mounted on an Unmanned Aerial Vehicle (UAV), and with a digital camera on the ground, combined with field measurements of herbaceous aboveground biomass.
The different measurements (UAV, camera, and HAB) were carried out on plots of 1 ha. The UAV images covered all the plot. Within the plot, subplots of 1m by 1m were laid out. The camera plot and the HAB measurements were carried out on these subplots ( Figure 1).
The positions of the subplots were marked out with a wooden triangle painted in white, or with a plastic bag on the ground.
For the different datasets, we first present the position of the plot and the subplots and the dates of the different measurements.
Second, we detail the protocol for the image acquisitions and image analysis (Structure from Motion).
Lastly, we present the herbaceous biomass field measurements.

| Design of the measurement setup
Data were collected on three spatio-temporal scales for comparison: 1. Temporal scale: images and HAB data were collected at the same locations over the growing season (corresponding to the wet season between June and October) (Figure 2a).
2. Landscape scale: measurements were done within a 6,800-ha area of enclosed savannah in northern Senegal (Figure 2a).
This design produced six datasets. Three datasets for the analysis of UAV imagery: Temporal UAV (TU), Landscape UAV (LU), and National UAV (NU); and three datasets for the digital camera on the ground: Temporal camera (TC), Landscape camera (LC), and National camera (NC).

| Temporal scale measurements
The temporal scale measurements were carried out at the Centre de Recherches Zootechniques in Dahra (Figure 1a), a research station in northern Senegal (15°21' N, 15°28' W) in the silvopastoral zone of the country referred to as the Ferlo region. The station covers an area of 6800 ha managed by ISRA (Institut Sénégalais de Recherches Agricoles). Mean annual rainfall over the last 50 years was 371.6 mm and the soils are mainly sandy (Ndiaye et al., 2014(Ndiaye et al., , 2015. Images were taken from 1 m 2 subplots every 10 days with the digital camera on the ground during the 2019 rainy season (end of July to October) and during the 2020 rainy season (July to October). In 2019, images were collected close to the Dahra research center office, and in 2020 they were collected around the Dahra flux tower site (Tagesson, Fensholt, Cropley, et al., 2015).
Both years, the images were collected along two perpendicular 1 km transects . Every 10 days, measurements were taken 1 m to the north of the last measurements (10 days earlier).

F I G U R E 1
Schematic representation of the plot and the subplot(s). The positions of HAB are the black squares measuring 1 m by 1 m. For the Unmanned Aerial Vehicle (UAV), the images were taken over the entire plot. The camera images were taken over the subplots. There were from 1 to 10 subplots per plot The UAV images were taken every two days during the 2020 rainy season, covering a one ha area. For the HAB measurements, three subplots of 1 m 2 were sampled every time in relation to the distance from the center of a tree: one under the crown of the tree, one at the perimeter of the tree crown and one outside tree influence (at a greater distance than the height of the tree) to be compared with HAB measurements collected on the ground from destructive sampling (see description below). The position of the plot changed each time. In total, the measurements were taken around six different trees of two different species (Balanites aegyptiaca and Vachellia tortilis).

| Landscape scale measurement
For the landscape scale measurements, data were collected from 38 plots of 1-ha close to the Dahra research station ( Figure 2a, the red squares). These plots had already been used in several studies of vegetation dynamics and were selected to be representative of the diversity of vegetation types within the research station (Ndiaye et al., 2014(Ndiaye et al., , 2015Raynal, 1964). Data were collected at the end of the rainy season (October 2018) when HAB was at its maximum. Images covering the full 1-ha plots were collected with the camera mounted on the UAV. For HAB, 10 subplots of 1-m² were measured along a 20 m line. The images with the digital camera on the ground were only taken for the first 1 m² subplot among the 10.

| National scale measurements
The data collected on a national scale were taken from 47 plots  For temporal and national scale data collection, the UAV was a Parrot Anafi with a Pix4D capture application (drone flight planning app) using a double-gridded flight plan that covered at least 1 ha. For data collection on the temporal scale, the height of the flight was 60 m with an overlap of 80% at low speed. All the flights were made early in the morning. Images were acquired in autofocus mode with an 80° angle.
For data collection on a national scale, we used the same doublegridded protocol proposed in the Pix4D capture application, but with a flight altitude of 80 m. Furthermore, the flights took place during the day at some point between 8 am and 7 pm.

| Ground-based camera
For data collection, a canon Ixus 180 (20 MP) was used for the landscape scale, while a Campark 4K Ultra HD 20 MP was considered for the national and temporal scales. For the landscape scale, we made videos using a digital camera above the 1 m² subplots. For the national and temporal scales, the cameras were moved above the grass along five 1-m transects. The camera had a 90° orientation in the direction of the ground. Close to each transect, we had an item with a reference height of 20 cm. For all the data collected, the videos were made at a height of 1 m from the ground.

| Structure from motion analysis
All the images were analyzed with Pix4D software (photogrammetry software) ( Figure 3). The videos from the digital camera on the ground were split into images. The split was such that there were between 200 and 300 images per subplot. We used the "3D maps" mode of the software (Figure3a and b), which produces an orthomosaic with a digital number for the red, green, and blue channels  For the outputs based on data from the digital camera on the ground, we extracted the data from the height reference and corrected the DSM of the different plots. From the CHM, we used the average height of all the pixels within the subplot (H mean ) and the maximum height of the subplot (H max ).
We used the digital numbers when calculating the indices. For the NU dataset, we also computed reflectance using the option in Pix4D software. The indices calculated for reflectance were closely correlated F I G U R E 3 Output from Structure from Motive (SfM) (a) 3D point cloud obtained from the Unmanned Aerial Vehicle (UAV) images, (b) 3D point cloud obtained from the camera images. (c) Orthomosaic obtained from the UAV images, (d) Orthomosaic obtained from the camera images. (e) Digital surface model obtained from the UAV images (The height is in meter above sea level), (f) Digital surface model obtained from the camera images (the height is in cm above ground) with the indices based on the digital number and we thereby concluded that the digital values could be used directly (Appendix S1).

| Collection of herbaceous aboveground biomass
Herbaceous aboveground biomass (HAB) was collected for each of the 1-m 2 subplots. The herbaceous vegetation was cut at soil level, as by Levang & Grouzis (1980). The fresh mass was weighed directly on site. Two different methods were used to estimate dry mass. For the temporal dataset, the whole-fresh sample was dried at 56°C and the mass was weighed daily. When the samples stopped decreasing in weight, we considered that weight as the dry mass. For the HAB data collected on the landscape and national scales, the large number of samples made it impossible to dry all the samples. We took a composite sample from each site with a known fresh weight, which we dried to obtain the dry matter content. We then applied the dry matter content found for the samples of fresh mass to obtain the dry mass of each sample (Diouf et al., 2015;Ndiaye et al., 2015).
The number of field measurement and the statistic of the fresh and dry mass across the different dataset is in Table 2.

| Predicting biomass from SfM outputs
For each of the six datasets even for the datasets with a low number of measurements (TC 29 and LC 35), we randomly separated the data into two groups; 2/3 were used for training (calibration dataset) and 1/3 were used for validation (validation dataset). The model calibration was done using a random forest algorithm with the field samples of fresh and dry HAB as response variables and SfM output as explanatory variables. Due to the unbalanced distribution of HAB (more low than high values), we used a root mean square transformation. The predictor variables were the mean and maximum heights (H mean and H max ), the red, green, and blue digital numbers, and the six vegetation indices. Each random forest used 500 regression trees, each time only including six randomly selected predictor variables from one-third of the calibration dataset. The percentage of explained variance was obtained by crossvalidation against the remaining 33% of the calibration dataset within the random forest. The importance of each predictor variable was estimated by taking the difference in prediction quality based on the mean square error when the randomly selected variable was included and not included in the regression trees. This difference was then averaged over all trees and normalized by its standard deviation.
The calibrated random forest models were used to predict fresh and dry HAB for the validation sites and compared these predictions to the field-observed HAB of the validation dataset. The agreements between model predictions and field observations were assessed using the root mean squared error (RMSE) and the relative RMSE obtained by dividing the RMSE by the mean of the measured values (RMSE R ). We also calculated the root median squared error (RMdSE) and the relative root median squared error (RMdSE R ).

| Comparing Camera and UAV output
The goal here was to compare Camera and UAV predictions on the same set of data for subplots where both UAV and camera images were taken. This was the case for 35 subplots of the landscape

| Combining UAV datasets
The random forest model was run on the combination of the three UAV calibration datasets (i.e., TU, LU, and NU). We obtained a percentage of explained variance of 78.65% for fresh mass and 68.09% for dry mass ( Table 2). The most important variables were the same for both fresh and dry mass (i.e., Gli, Vari, Exg).

| Combining camera datasets
For the model based on all the camera data, the percentage of explained variance for fresh mass was 47.81%, and 35.58% for dry mass. The three most important variables were Gli, Exg, and H mean for fresh mass, and H mean , Gli, and Vari for dry mass.
The R² for fresh mass when combining the camera datasets

| Comparison of camera and UAV data
For the variance partitioning of fresh mass (Figure 6a), 58% of the variance was explained by both types of data, 14% only by the camera and 2% only by the UAV. The residuals (unexplained variances) amounted to 26%. For variance partitioning for dry mass with 33% of residuals (Figure 6b), 50% of the variance was explained by both types of data, 12% only by the camera and 5% only by the UAV.

| Assessing herbaceous aboveground biomass using SfM outputs based on RGB
Our results showed that Structure from Motion (SfM) data obtained from RGB images can be used to assess HAB of the herbaceous layer with relatively good accuracy. For fresh mass, we had errors of around 150 g m − ² across the different datasets (around 25% of relative error), whereas for dry mass, the errors were around 60 g m − ².
However, there is some room left for improvement to be discussed below. SfM outputs are from two types of information: color indices ( Table 1) and height variables derived from the DSM and DTM.
Regarding the most important variables, for the color indices the color index based on all three colors (Gli, Exg, Vari indices) proved especially useful. At least one of these variables proved to be among the three most important variables in the 24 models we tested. The green leaf index (Gli) was amongst the three most important variables in 13 models tested and the most important variable in eight models. This is supported by the literature, where the same results were found (Bendig et al., 2014;Lussem et al., 2018;Possoch et al., 2016). For most of the datasets produced, the images were taken at different times of the day and therefore with different light conditions. A color grid was used as part of four of the six datasets, but the correction made with these color references proved to induce more errors at the end (results not shown) and therefore this preprocessing step was omitted from the results presented here. Furthermore, the index based on reflectance data was found to be less related to field measured data than indices based on digital numbers. This result shows that the color indices were related to biomass. Here, we only based our work on RGB images of course, but many studies have used near infrared. (Candiago et al., 2015;Wijesingha et al., 2020). Indeed, it is well known that Infrared-based indices are linked with the vegetation (Rouse et al., 1974). Increasing the spectral resolution of UAV captors could also be a solution to increase the prediction. However, the cost of a NIR UAV or camera is high, but it would be interesting to compare RGB indices and NIR indices.

The height variables (maximum and mean height of the subplot)
obtained from the 3D model were within the top three variables in nine out of 18 models, especially for the dry matter model (six out nine models). The height variables were more important in the models based on camera outputs. Spatial resolution was found to be quite different between the two tools, whereas for the UAV the ground surface distance (GSD) was around 3 cm (a few mm for the camera).
The height of the grass in the Sahel region is generally around 40 to 50 cm at maximum vegetation growth (Boudet, 1984). GSD may not be precise enough for the height of grass and its variation within the natural vegetation. In other studies, using UAVs on grasslands, the variables from the 3D model were found to be more predictive (Lussem et al., 2018(Lussem et al., , 2019. In some studies, the measurements were taken on management experiments that created high variability in grass heights between the different plots. Their flight height was lower than in this study. We chose a height of around 80 m due to the presence of trees. It was not that the height of the tree could directly impede a flight altitude of 20 m (the highest tree was around 15 m tall), but because the 3D model of the trees would have been sub-optimal and may have impeded the creation of the mosaics.
Furthermore, we also used the Landscape and National scale images to work on tree communities . For the camera, an area of bare soil was close to the subplot. Indeed, we cut the grass close to the subplot. This bare soil created a difference in height variability and the spatial resolution was much better. For Errors could also come from the manual measurements of the herbaceous vegetation (destructive sampling). For example, the percentage of variance explained was generally higher for fresh mass than for dry mass, except for the analysis covering the temporal scale. One reason might be that in the datasets of the temporal scale analysis the dry mass was measured for all samples. For the other four datasets, the dry matter content was only measured for a composite sample. The composite sample may have created some bias in the measurement (Diouf et al., 2015).
One problem is also the positioning of the measurements, especially for the UAV. We used markers on the soil (wooden triangle or plastic bag), but some inaccuracies still remained. Ideally, it would be better to have a 1 m² subplot with a visible border.

| Comparing camera and UAV tools
We investigated using Structure from Motion (SfM) from two different types of sensor systems: Unmanned Aerial Vehicle (UAV) and a ground digital camera. Overall, we found higher percentages of explained variance for the UAV than for the camera. However, when we compared the same set of subplots for both tools ( Figure 5), the camera outputs were slightly better than those of the UAV. We surmise that the difference observed between the two tools was due to a better assessment of vegetation height by the camera.
Although both tools are based on the same concept, they have quite different uses. UAVs are now frequently used for vegetation monitoring. Our work confirmed that UAVs, especially when equipped with a low-cost RGB sensor, can be used to adequately quantify the phytomass of herbaceous layers in a savannah ecosystem. UAVs can be used to map areas from 1 ha up to 20-30 ha (potential area with one UAV battery) depending on the flight plan and UAV characteristics. The spatial scale of this approach is well suited to experiments on permanent grasslands (Lussem et al., 2019). For many temperate ranching systems, UAVs could also be used to assess biomass on a plot scale to evaluate the available livestock feed.
However, in Sahelian countries, as in many pastoralism areas, rangelands are used as a common resource exploited by several farmers across large areas. UAVs could be used as an intermediate tool to close the spatial gap between field observations and satellite images. Many UAV studies have shown that woody vegetation can also be characterized using UAVs Mayr et al., 2018). Trees are an important part of the savannah ecosystem, even with their low density. One interesting feature of the UAV tool is the ability to evaluate both strata of vegetation . It can also be used to assess the spatial struc-  (Haula & Agbozo, 2020).
Indeed, the use of UAV technology is limited by law. In many countries, some areas are restricted, such as populated areas, airports, or other sensitive areas (e.g., military) and for some countries the use of UAVs is totally forbidden and/or detailed legislation has not effectively been implemented. Furthermore, UAV licenses are required in many countries, which means that UAVs cannot be used by just anyone, and that potential users will have to take training courses to obtain the license.
Using SfM approaches with a simple digital camera has been much less studied than with UAVs (Andújar et al., 2018;Cooper et al., 2017), even though it is easier and cheaper to use a camera. In the study by Cooper et al. (2017), SfM was obtained from photography and here we used a video to simplify the field work and to be able to subsequently select the number of images implemented in PIX4D software.
The digital camera has advantages for use on a fine spatial scale, due to its higher definition (a few mm as opposed to 3 cm) as compared to UAV systems. In cases of low phytomass, such as early in the growing season, or for very low and sparse herbaceous layers, individual grasses can be easily identified (Figure 1). Structure from Motion has already been used on individual plants Andújar et al., 2018)  be assessed, such as the phytomass or the mean height of the layer. (Bossoukpe et al., 2020).
Another interesting feature of the camera approach is its ease of use, and the tool can be used to develop a herbaceous growth monitoring programs across a large scale (national scale). In Ireland, a participatory monitoring program has been in operation for several years, where farmers across the country estimate grass growth (Hanrahan et al., 2017). The data are available online and can be used by farmers or the authorities. In Sahelian countries, this kind of data could greatly enhance the current monitoring framework (Diouf et al., 2015).
However, compared to UAVs, using a digital camera has several drawbacks. Of course, a digital camera cannot be used to produce maps and the orthomosaic is not geotagged. Moreover, a height reference is needed. For the UAV, several subplots can be taken for a single flight, whereas for the camera each sample requires an SfM process and the number of pix4D projects needed for calibration will be larger (one per sample).

| Using the SfM methodology across spatiotemporal scales
Our results showed that Structure from Motion (SfM) can be used on the three different scales tested here. Regarding the different scales and datasets, the range of the biomass masses was quite similar between the different datasets ( Table 2). Between the six datasets, the percentage of variance was clearly higher for the datasets including more measurements. All in all, we observed similar errors across the three different scales and showed that SfM can be used in all instances of scales.
On the temporal scale, we tested SfM use during the growing season. The growing season in the Sahel region is very short and erratic during the few months of the wet season and UAVs and cameras can be used to monitor grass growth with a high measurement frequency. On an experimental design with differences in management or water input on small plots, such as (Cisse, 1986;Hiernaux & Hérault, 2003;Hiernaux et al., 1995;Hiernaux & Turner, 1996), SfM could be used to monitor the dynamics of the herbaceous layers indirectly with high frequency.
Here, we only tested the method during the wet season and on green grass. However, dry grass during the nine months of the dry season is an important source of animal feed in the region. The continuous monitoring of this remaining dry biomass is key for livestock management in pastoral systems. Remote sensing work based on the soil tillage index using shortwave infrared information has been carried out to evaluate dry biomass . Our analysis of monitoring on a national scale involved a dataset collected in the northern and southeastern parts of Senegal. In addition to selecting plots based on topography and management type, the plots were located along a rainfall gradient (from an average of 200 mm per year to more than 1000 mm per year), and a soil type gradient from sandy soils to ferralitic soils was also considered.
A national scale monitoring program is important for livestock policy, and indeed in several Sahelian countries the available biomass (herbaceous layers and foliage of woody species) is monitored on a national scale taking a sample-based approach. The aim is to evaluate the national fodder balance and eventually develop local and national policies (pastoral mobility, distribution of animal feeds, etc.) (Touré & Ickowicz, 2014). Our results showed that SfM tools could be used across Senegal and possibly other Sahelian countries. For the UAV approach, one limitation is the southernmost region where the tree cover can be dense, and only the grass biomass outside the tree cover can be monitored.
UAVs or other very high resolution (VHR) remote sensing systems could be used to simplify scaling-up from field measurements and conventional remote sensing systems (Taugourdeau et al., 2014).
Indeed, a hybrid model combining fewer measurements could replace actual field measurement programs that require an extensive protocol to assess the heterogeneity of the biomass layer, as already proposed for other communities (Kattenborn et al., 2019).
The camera approach could also be used on a national scale. One interesting feature of a camera monitoring protocol is that it requires few qualifications and equipment and one way forward will be to use this tool as part of a participatory observatory. Different observers across the country could take RGB images using a digital camera/ smartphone, which would be analyzed to monitor biomass across the country.

| Combining the use of datasets
We obtained better results for biomass estimation when more sampling data were used. One way forward would be to combine several datasets to build a large calibration dataset. However, our results showed that prediction was not better and generally performed slightly worse when our calibration was based on a combination of several datasets, as compared to when calibration was based on a single dataset (Figures 3 and 4). This result was found even for the dataset with few observations, such as the "Temporal Camera" (29 plots) and the "Landscape Camera" (35 plot) approaches. For these two datasets, the random forest and the separation in calibration set and validation set are necessarily adapted due to the low number of measures but we used the same methodology for all datasets.
For both, the percentage of variance of the random forest was quite low, but the prediction on the validation dataset was not markedly different from the other datasets. The difference between the combined datasets and single datasets was less pronounced when the combination of datasets was made between UAV and camera, separately. The small difference between the results based on calibration within a dataset and calibration made by combining the datasets from the same tool showed that tool-specific calibrations were quite generic. It can be hypothesized that the calibrations developed in this work could be used directly on new datasets of each type.
We found small differences between the makes of UAVs (Spark or Anafi) and cameras used (Campspark and Canon) and on the image acquisition protocols (altitude, type of flight plan).
It would be interesting to produce more data using UAVs on herbaceous plants involving different types of vegetation and protocols. The goal of such work would be to further test the robustness and generalness of UAV biomass relations for African savannahs.

| CON CLUS ION
We showed that RGB images used for structure from Motion (SFM) can be used to assess aboveground biomass of the herbaceous layers in a Sahelian savannah ecosystem. We tested SfM processing for biomass estimates using data from Unmanned Aerial Vehicles (UAV) and ground digital cameras and on three scales. The three different scales included spatial scales within a landscape or across a country and, lastly, the temporal scale measuring continuously over the growing season. Our results showed that this approach could be adequately used on all the different scales and testing of the methods also took into consideration variations in land management and soil types for a rainfall gradient from arid to sub-humid conditions. The ground digital camera data were found to be slightly more accurate than the UAV data, but UAV measurements have the advantage of facilitating the scaling-up of field measurements to areas of several ha. Herbaceous aboveground biomass mapping by SfM can be used to undertake studies of ecological processes at a high spatial resolution and provides an image source building an important bridge in spatial scales between field observations and satellite remote sensing images.